The Radiative Kicked Oscillator: A Stochastic Web or Chaotic Attractor ? 
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A relativistic charged particle moving in a uniform magnetic field and kicked by an electric field 
is considered. Under the assumption of small magnetic field, an iterative map is developed. We 
consider both the case in which no radiation is assumed and the radiative case, using the Lorentz- 
Dirac equation to describe the motion. Comparison between the non-radiative case and the radiative 
case shows that in both cases one can observe a stochastic web structure for weak magnetic fields, 
and, although there are global differences in the result of the map, that both cases are qualitatively 
similar in their small scale behavior. We also develop an iterative map for strong magnetic fields. In 
that case the web structure no longer exists; it is replaced by a rich chaotic behavior. It is shown that 
the particle does not diffuse to infinite energy; it is limited by the boundaries of an attractor (the 
boundaries are generally much smaller than light velocity). Bifurcation occurs, converging rapidly 
to Feigenbaum's universal constant. The chaotic behavior appears to be robust. For intermediate 
magnetic fields, it is more difficult to observe the web structure, and the influence of the unstable 
fixed point is weaker. 
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I. INTRODUCTION 

Zaslavskii et al |jj studied the behavior of particles in 
the wave packet of an electric field in the presence of 
a static magnetic field. For a broad wave packet with 
sufficiently uniform spectrum, one may show that the 
problem can be stated in terms of an electrically kicked 
harmonic oscillator. For rational ratios between the fre- 
quency of the kicking field and the Larmor frequency 
associated with the magnetic field the phase space of the 
system is covered by a mesh of finite thickness; inside 
the filaments of the mesh, the dynamics of the particle 
is stochastic and outside (in the cells of stability), the 
dynamics is regular. This structure is called a stochas- 
tic web. It was found that this pattern covers the entire 
phase plane, permitting the particle to diffuse arbitrarily 
far into the region of high energies (a process analogous 
to Arnol'd diffusion (|]). 

Since the stochastic web leads to unbounded energies, 
several authors have considered the corresponding rela- 
tivistic problem (for a work which can be related to the 
relativistic stochastic web see ||). Longcope and Sudan 
j| studied this system (in effectively 1^ dimensions) and 
found that for initial conditions close to the origin of the 
phase space there is a stochastic web, which is bounded 
in energy, of a form quite similar, in the neighborhood 
of the origin, to the non-relativistic case treated by Za- 
slavskii et al. Karimabadi and Angelopoulos J5|] studied 
the case of an obliquely propagating wave, and showed 
that under certain conditions, particles can be acceler- 
ated to unlimited energy through an Arnol'd diffusion in 
two dimensions. Since an accelerated charged particle ra- 
diates, it is important to study the radiative corrections 
to this motion. We shall use the Lorentz-Dirac equation 
to compute this effect. 

We compute solutions to this equation for the case of 



the kicked oscillator. Under the restriction of weak mag- 
netic field, at low velocities, the stochastic web found 
by Zaslavskii et al JIJ occurs; the system diffuses in 
the stochastic region to unbounded energy, as found by 
Karimabadi and Angelopoulos ||. The velocity of the 
particle is light speed limited by the dynamical equa- 
tions, in particular, by the suppression of the action of 
the electric field at velocities approaching the velocity of 
light §. 

For the case of a strong magnetic field case (a case 
which can not be examined in laboratory) the stochastic 
web, which may occur in the weak magnetic field case, 
does not appear and is replaced by a rich chaotic behav- 
ior. In this regime, the particle does not accelerate to 
infinite energy; it limited by the boundary of the chaotic 
attractor. 



II. MODEL 

In the present study we will consider a charged par- 
ticle moving in a uniform magnetic field, and kicked by 
an electric field. The effect of relativity, as well as the 
radiation of the particle, will be considered. We restrict 
ourselves to the "on mass shell" constraint (keeping con- 
stant x^x,. = 



2 )Q 



The fundamental equation that we use to study radia- 
tion is the Lorentz-Dirac equation jsj , 

e I 
m i M = -XyF^ + 7 m (i' M - — x^x v x v ), (I) 
c cr 



*The more general case which includes also the possibility 
"off mass shell" motion, is discussed in Ref. 11711. 
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where 70 — - — 



. — 6.26 x 10~ 24 sec. The indices and v 
indicate the coordinates t, x, y, and z (or 0, 1, 2, 3), and 
the derivative is with respect to r, which can be regarded 
as proper time in the "on-mass-shell" case. F^ v is the an- 
tisymmetric electromagnetic tensor. The first term of the 
right hand side of Eq. (Q) is the relativistic Lorentz force, 
and the second term is the radiation-reaction term. Note 
that the small size of the radiation coefficient (70 <C 1) 
leads to a singular equation that requires special mathe- 
matical treatment, as well as physical restrictions, as will 
be discussed in the succeeding sections. 

Following Zaslavskii et al Q , the magnetic field is cho- 
sen to be uniform in the z direction, and the kicking 
electric field is a function of x in the x direction, 

B = (0,0, B) 

OO 

E(x,t) = (f(x) 5(t-nT),0,0). (2) 

n— — 00 

Originally, Zaslavskii et al chose a uniform broad band 
electric field wave packet which can be expanded as an 
infinite sum of (kicking) (^-functions (u>o is the frequency 
of the central harmonic of the wave packet, ko is the wave 
number of the central harmonic, and Auj is the frequency 
between the harmonics of the wave packet) 



E x = E(x, t) = -E 



E 



siri(kox — u>ot — nAujt), (3) 



which, for luq = 0, becomes 

00 

E(x,t) = -E Tsm(k x) ^ 5(t-nT), (4) 



where T = Eq. (g) is the generalization of Eq. 

(||), with an arbitrary function fix) instead of the sine 
function of Eq. (||). 

In order to introduce a map which connects between 
cycles of integration which start before a kick and end be- 
fore the next electric field kick, one has first to integrate 
over the 5 electric kick and then to integrate the equa- 
tions of motion between the kicks, where only a uniform 
magnetic field is present and there is no electric field up 
to the beginning of the next kick. 



III. THE MOTION OF A CHARGED PARTICLE 
IN A UNIFORM MAGNETIC FIELD 

A. An integration with respect to the proper time r 

In the case of a uniform magnetic field in Eq. (||), Eqs. 
(|l|) reduce to three coupled differential equations, 

ci — x — 7o(i'o ?ioR) 



X = X\ = —Q.X2 + 70(^1 TrXiR) 

y = x 2 = flx 1 + 70(2:2 - \x 2 R), 



(5) 



where — R — x\+x\~ Xq, and e = — e (charge 

of the electron). 

Using a complex coordinate u = x + iy, Eqs. (|5|) 
can be written as, 



1 



t = 7o(i - -ntR) 
u = iflii + jo(u ^iiR). 



(6) 



It is very convenient to use the hyperbolic coordinates, 
i = cosh q 

x = c sinh q cos <f> (7) 
y = c sinh q sin (f> 

since the constraint (see 0] for further discussion), 

• u ■ ■ 2 1 -2 2/2 2 /n\ 

x^x^ ~ x + y — c t = —c , (8) 

is then automatically satisfied. The complex coordinate 
u then becomes, 



= csinh qe 1 ^ 



(9) 



and the equations of motion Eqs. m) can be written as, 



q = — 7o^ 2 cosh q sinh q + joq 
(p = Q + 27o0gcoth<3 , + 7o0- 



(10) 
(11) 



As pointed out earlier, in the case of a singular equa- 
tion (such as Eqs. (10) and (pi])), one has to consider 
physical arguments as well as mathematical arguments. 
There are several suggested methods to avoid the "run 
away electron" problem ||. One of the most frequently 
used, especially useful for scattering problems, assumes 
JlOt that the particle loses all its energy after a sufficiently 
large time, and the run away parts of the solution are set 
to zero; the equations of motion can then be written as 
an integral equation. Another method, which permits 
the study of problems with only non-asymptotic states 
jnj, uses an iterative singular perturbation integration 
that leads to the stable solution. Since the integration 
we must use is over a finite time, it is impossible to use 
an asymptotic condition (e.g. |u| — > at t — > 00). The 
approach of Sokolov and Ternov H is more suitable for 
our purpose, since it can be implemented for bounded 
time integration, and the mathematical solution is sim- 
ple. In this method, in the first step, the perturbation 
terms in Eq^ ( p"l| ) are disregarded (as in the iterative 



scheme of |11 ), and then the resulting equation is ex- 
actly integrated. In the second step, the singular term 
on the right hand side of Eq. ([l(]) is not considered; us- 
ing </>(t) from Eq. ([ll]), the equation including just the 
first term is integrated. The solution is 



2 



<j) = flT 

(3{t) = tanhq = ^° e -7on 2 r = 
c 



(12) 



where /3q = 21 is the actual initial velocity divided by 
c, and To = 7o q^ is the decay time for the energy of the 
particle. 

To get an estimate for the radiation during one cy- 
cle, consider the maximal uniform magnetic field that 
can achieved today in a laboratory, which is[] O(10T) = 
O(10 5 gauss). If, for example, B = 10 5 gauss, then 
il = 1.76 x 10 12 — , and r = 5x 1CT 2 sec. Thus, it is clear 

sec 5 u 1 

from Eq. ( |l2| ) that the particle makes J^t ~ 10 10 cycles 
before it decays to = of its initial velocity. In other words, 
the energy loss during one cycle in very small, and since 
in our problem, the time T between the kicking is of the 
order of the period, the energy loss between consecutive 
kicks is very small. 

The time between the kicking is measured according 
to the observed time along the motion, At = T. Thus, 
it is necessary to find the corresponding At. It follows 
from Eq. (g) and Eq. @ that, 

i = cosh q = — (13) 

The solution of Eq. (H) can be obtained by a elementary 
integration, 



1 + VT 



-Ci 



t = tq In 



where a = foe T ° . After some algebra one gets, 



At 



/l ^ + q + y/T^V^ 

TOn U (l + v /T-7f) e T/ro 



(14) 



(15) 



In Fig. [l] we present the behavior of the function 
AT TT) ' ^ can ^ e seen ^ rom Fig- II * na * m an y case ; 



At < T: 



(16) 



this implies that the time difference according to the 
proper time r is always less then the time difference ac- 
cording to observed time t. This fact is corresponds to 
the well known relativistic "time-dilation" phenomenon. 
The inequality Eq. ( [l6|) can be also illustrated by con- 
sidering the two extreme cases, 



r«T 

T»T 



=> At w TJl 



Pi 



2 t 



0(T 3 ) 



At « T + t In 



i + ^r~w 



(17) 
(18) 



^It is possible to obtain a short-lived magnetic field of 
O(10 6 gauss) in the laboratory. 



For strong magnetic field (i.e. the case of Eq. (Jig)) the 
numerator is less than (or equal) to 2 and thus the In 
function will give a negative number; in this case the in- 
equality Eq. ([16|) is clearly achieved. For a low magnetic 
field Eq. Q behaves like a parabola; the slope at T = 
is \/\ — 01 which is less than one, and thus the inequal- 
ity Eq. (|16|) is satisfied. In the present study we will 
consider separately the cases of weak magnetic field (Eq. 
( |l7| ) and strong magnetic field (Eq. (|l8|)). 



Fig. 1 



Ax/x„ 




FIG. 1. The function At( ^ /to) versus T/r . Four typical 
cases are shown (/3 = 0, /3 = 0.7, /3 = 0.9, /3 = 1). All 
possible curve lines between the graphs of /3o = and /3q — 1. 
The dashed line indicates the asymptotic behavior of the case 

A) = l. 



B. An integration with respect to the observed time 



In the previous subsection we have shown the approxi- 
mate solution for a charged particle moving in a uniform 
magnetic field. However, as mentioned above, the kick- 
ing of the electric field from Eq. (^) is with respect to 
the observed time t, and it is therefore convenient to inte- 
grate the equations of motion (Eqs. ( |l0| ) and (p~T|) ) with 
respect to t. 

The arguments that were used in order to obtain Eqs. 
( |l2| ) , as well as the implementation of the chain rule on 
Eqs. © and @ by the use of Eq. (@), lead to the 
following equations of motion 



dq 
~dt 



■ sinh q 



dt cosh q 
A simple integration of Eq. ( |l9| ) yields, 



(19) 
(20) 



3 



tanh ■ 



e Cl e , 



where C\ < 0. Substitution of Eq. in Eq. 



(21) 



d0 _ 
~dt ~ 

and the solution is 



fitanh( Ci), 



to 



= £!Tolncosh( Ci) 



C 2 



(22) 



(23) 



where Cj = cosh (-gj-)j and C 2 = tan 1 (yo/x ) + 
r2roln/3o- Returning to the original x, y coordinates (us- 
ing Eqs. (0)), the equations of motion become, 



dx cos(f2r ln(cosh(i/r — Ci)) + C 2 ) 
dt cosh(t/r — Ci) 

dy sin(f2T ln(cosh(t/r — Ci)) + C 2 ) 
dt cosh(t/r — Ci) 



(24) 



It is possible to write Eq. (E 
by using the relation 

cosh(— — C\) = — cosh — 
to Po T 



in more convenient way, 



-Vl-/3 2 sinh- (25) 



Tf) 



Eqs. ( p4| ) then become 

(f)o cos «-(f) si 



<i?; 
dt 



sin a 



cosh -i- + \/l - /?n sinh — 

cos a 



(f ) sina 



( rf «)o 



cosh — + */l - sinh ^ ' 



(26) 



in a Ta ylor series and then to integrate. The expansion 
of Eq. @ is 



1 ^i 2 



1 



Z TQ O Tn Z Tq 



(28) 

The first order expansion (according to t/r ) of Eqs. j2^ ) 
is 



dx 
~dt 



dy 
dt 



1 



1 + II 

7 t 



cos(-t+-ft/3 2 -) 
i 7 2 T 

sin(-t + -fi/? 2 -) 

I 7 -s t _ 



7 to 



dx\ , ,n l^^i 2 ^ 



where 7 = 1 
approximation 



Using Eq 



cos(—t - 

o 7 

l t 



1 t 2 

2 Tq 



(29) 



1 — Pq. Since -C 1, one can use the 



1+1± 

7 to 



1 t 

7 To' 



(30) 




the actual velocities, dx/dt and dy/dt, 
from Eq. (I29J) can be integrated and expressed by elemen- 
tary Fresnel functions. However, the effect of radiation 
is due to the term (which multiplies the sine and 

cosine functions), and the ^Q/3q^ term is not essential 
since the major offset from the unstable fixed point is 
due to the —t term. 

7 

Under the above assumptions the resulting equations 



where 



£Itq ln(cosh — (1 
to 



/Sgtanh— )). 
to 



(27) 



It is necessary to integrate Eq. (|26j) since the value of 
x is used in the electric field kicking in Eq. (|J). There 
is no analytical solution to Eq. (26), and a numerical 



integration must, in general, be performed. In the suc- 
ceeding sections we will study separately the case of weak 
magnetic field (i.e. T/tq <C 1) and the case of strong 
magnetic field (i.e. T/r ^> 1). 



IV. DERIVATION OF THE MAP - WEAK 
MAGNETIC FIELD 

A. The approximated solution for the weak 
magnetic field case 

When dealing with weak magnetic field we restrict our- 
selves to T/tq <C 1, and it is possible to expand ^ and 4| 



=■) =11- 



T 

7 T o 



dx\ ,n . 

— cos-T 
dt / n 7 



^ = 1- 





T 

7 r o 



dx \ . n . 

— sin-T 
at Jo 7 



One therefore obtains 

dy 
~di 



sin(-T) 
o 7 . 



cos(-T) 
o 7 . 



7 



« 2 t 



daA 
dt /, 



dx\ 
dt J, 



a; . 



(31) 



(32) 



(33) 



The exponential decay from Eqs. ( |24| ) is replaced by 
linear decay; by squaring Eq. (|3|) and Eq. (^2|) one 

obtains the linear relation v(T) = 



T/it ). 



.,2 _ 



«b(i 



4 



B. Integration over a S electric kick 

As pointed out in the previous section, the Lorentz- 
Dirac equation (Eq. ([[])) is a singular equation. The 
electric field which is used in this paper was expanded to 
a sum of S functions (Eq. (0)). In that case, Eq. (|l|) can 
not be integrated over the kick by regular treatment, and 
some approximation for the S function should be consid- 
ered instead. However, since the radiation has its effect 
in the direction of the relevant coordinate (in our case it 
is the x 1 — x coordinate), and because of the infinites- 
imal time interval, it is possible to replace the kicking 
parameter, K, by an effective kicking strength, which 
should be a smaller number than the original one. Thus, 
the nature of the solution for the system will not change 
due to the radiation during the kicking. The constant 
magnetic field during the kick is not considered since the 
integration is over an infinitesimal time interval. 

Under the above assumptions, in the neighborhood of 
the kick, the Lorentz-Dirac equation, Eq. (Q), can be 
written as, 



d ( dx 
dt \aW 

d_(dy_ 
dt \ dr 



= f(x) £ 5{t - nT) 



71— — OO 



= 0. 



(34) 



Integration over the S function yields, 



dr _ 



(35) 



where the + sign indicates after the kick, and the - sign, 
before the kick. Following Ref. Q, we choose, 



fix) 



eET 



■ sin kx = K sin kx. 



m 



(36) 



Using the fact that 4- — 7^ one obtains the following 
relations, 



7+ = 



\ 



1 

! + - 

c z 



17 +/(*«) 



. 7 2 (dy 
J - 1 dt 



dt J 



1_ 

7+ 
1_ 

7+ 



dt J 



+ —f(x n ) 

7+ 



(37) 



n+l V 7+^0 



n+l 



dx\ fdy\ , n . 

z cos( — T) + ( -r sm — T) 
dt) + V ; \dt) + V 7 + ' 

T 

1 I 

7+ 



dx\ , n , fd y \ , n s 

- sin(— T)+ (4) cos(-T) 



(38) 



dt J 



7+ 



dt 



7+ 



(39) 



~ 7+ 



7+ 



Wt 



n+l 

dx\ f dx 



dt 



n+l 



dt 



(40) 



Notice that the minus signs in Eq. fl37|) refer to the 



,"1 



points of the map. 

It is possible to return to the non-radiative limit by 
letting ^- — > 0. In that case the map is equivalent to the 
map which was derived by Longcope and Sudan Q . The 
nonrelativistic limit, which was derived by Zaslavskii et 
al pi, is achieved by letting — > and c — > 00. 



C. Analysis 

Up to the derivation of Eqs. (|38|)-(^0|) we have passed 
three stages, (a) the velocities were calculated using an 
exact solution and the position x can be calculated by 
numerical integration (Eqs. (|2^)), (b) the velocities were 
approximated by Eqs. ( p9[ ) and the position x by the 
elementary Fresnel function, and (c) the velocities were 
approximated by Eqs. (|3^)-(39) and the position x (Eq. 
(|40|)) was derived using elementary integration. Among 
all the combinations of the solutions for the velocities 
and the position x we choose the combination of the ex- 
act solution of velocities (case (a), Eqs. (p6|)) and the 
exact solution for the position x (stage (c), Eq. ([40|)). 
The value of x was selected from the third stage of our 
derivation since it enters the equation just in the kicking 
term, and there it affects only slightly the phase. Thus, 
one does not expect that this fact can change the typical 
behavior of the particle. The above analysis is summa- 
rized in Fig. ^. 



Returning to the initial charge e = — eo and thus Q — > 
— J7, the map which connects the velocities from just be- 
fore the kick to the next kick (by the use of Eqs. (|3~i[)- 
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Fig. 2 



(otherwise, the particle will not diffuse, and will not cre- 
ate a web structure), 



exact solution 



semi-exact 
solution 



approximated 
solution 




numerical 
integration 

by elementary 
Fresnel functions 

exact 
solution 



FIG. 2. A diagram which present the all possibilities for 
the map construction. The wider line indicates the chosen 
combination. 



The radiative map is then 



dx\ 



.dt) + 



(*). 



sin a 



n+l 



n+l 



X n +1 



cosh ^ 

_ ( dx} 
\ dt / 



1 — (3 + sinh 



cos a 



cosh — 

TO 



7+ 

n 

7+ 
Wt 



01 sinh ± 



(41) 



(42) 



n+l 



where the values immediately after the kick are as in 
© and a is defined in Eq. © (a — > — a since e = 
. Note that for our limit (T/tq -C 1) there is not a { 
difference between the map given by Eqs. (^8|)-(^0|) 
Eqs. 



(43) 

Eqs. 

-eo) 
;reat 
and 



V. RESULTS - WEAK MAGNETIC FIELD 

In order to obtain a web structure, there are two con- 
ditions that have to be fulfilled. In the nonrelativistic 
case, the first condition is that the ratio between the gy- 
ration frequency, Q, and the kicking time, T, is a rational 
number, a condition which can be expressed as follows, 



QT = 2tt 



P 



dx\ 

~dt)o 

dy 
dt 







(2n + l) 



(45) 



where n is an integer number. The symmetry of the 
web is determined by p and q. If, for example p = 1 
and q = 4, the particle is kicked four times during one 
cycle, and thus, the symmetry of the web will be a four 
symmetry Jl|. 

In the relativistic case [Q, the above conditions 
ar e sligh tly different, because of the additional factor 
yl — /3q which multiplies the OT term. However, if the 
initial velocities are small, i.e. vq <C c, the additional 
factor is close to 1, and thus the structure of the web 
should not change. In the radiative and relativistic case, 
as well as the non-radiative relativistic case, conditions 
and (Ej|) become, 



= =0 



dx 

Tit 



(46) 



For sufficiently large initial velocities the above condi- 
tions do not hold anymore, and the web structure is not 
observed. 

In this section we will compare (qualitatively; a quan- 
titative treatment for the diffusion rate will be studied 
elsewhere) the diffusion of the non-radiative particle and 
the radiative particle. Intuitively, one would expect that 
a radiative particle will diffuse more slowly than a non- 
radiative one, since the radiation effects act like friction, 
and are thus expected to "stop" the particle. However, 
this naive expectation is not true, as we will demonstrate 
below. 

In order to investigate the above assumption, we have 
chosen a four symmetry structure (q — 4). We used 
the same initial conditions for all cases; just the kicking 
strength K has been changed. The parameters values 
and the initial conditions were B = 10, v x _o — 0, v Vt o = 
4 x 1CP 5 , k = 10, and the number of iterations was N = 
10 6 . In Fig. ||, we present the results as follows: the 
first column is the non-radiative case, the second column 
is the radiative case, and the third column is the total 
velocity versus the iteration number (the circles indicates 
the non-radiative case and the squares the radiative one). 
In the first row K = 1 x 10~ 5 , in the second K = 2 x 10~ 5 , 
and in the third K = 4 x 10~ 5 . 



(44) 



where p, and g, are integer numbers. Secondly, one must 
start from the neighborhood of the unstable fixed point 
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Fig. 3 
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FIG. 3. Different types of diffusion behavior as described 
in the context. 

As can be seen clearly from Fig. |[ the web structure 
is valid for the radiative case (top panel, second column) . 
The web structure is also exists in the other cases; the 
web cells are very small and they are difficult to see. Al- 
though the radiative particle and the non-radiative par- 
ticle start from the same initial conditions they behave 
differently. 

In some cases the diffusion rate of the radiative case is 
larger than the diffusion rate of the non-radiative case, for 
example, in the top and the bottom panels after 9 x 10 5 
iterations the radiative particle reaches a higher velocity 
than the non-radiative one. On the other hand, in the 
middle panel the radiative particle is slower then non- 
radiative one. Obviously, it is impossible to draw any 
general conclusion for the question of whether a non- 
radiative particle is faster then the radiative one. For 
this a more systematic approach should be carried out. 
This would be beyond the scope of the present paper and 
will be considered elsewhere. 

When the velocity of the particle reaches close to the 
velocity of light it actually stops growing (according to 
evolution in the time t) although it increases its energy. 
In Fig. |l| we show an example for that behavior. The 
parameters values are : B = 10, v x , — 0, v y fi = lx 10~ 5 , 
k = 10, K = 0.01, and the number of iterations was N — 
1 x 10 6 (every 10 th iteration was plotted). Fig. ^a and 
Fig. ^c are the non-radiative cases, while Fig. Bp and Fig. 
^d are the radiative cases; in Fig. ^a and Fig. Ib the ma' 
of v x versus v y is presented, while in Fig. ^c and Fig. 
the total velocity v versus the iteration number is plotted 
(every 10 th iteration was plotted). As seen from Fig. ||, 
the particle accelerates quickly to high velocity, and most 
of the time stays with this high velocity, approaching 
light velocity asymptotically, while increasing its energy 
to infinity. 




2x10° 6x10° 10x10" "'" 2xl(f 6x10° 10x10 

n n 

FIG. 4. The behavior of a particle with high velocity, (a) 

the non-radiative map, (b) the radiative map, (c) the total 
velocity of the non-radiative case, and (d) the total velocity 
of the radiative case. 



VI. DERIVATION OF THE MAP - STRONG 
MAGNETIC FIELD 

A. Introductory remarks 

In the preceding section we discussed the realistic case 
of weak magnetic field. In this section we will consider 
the other limit, i.e., the strong magnetic field limit. De- 
spite of the difficulties in producing a strong, uniform, 
and stable magnetic field, one can assume that such fields 
do exist near astronomical objects, such as, neutron stars. 

When increasing the magnetic field strength the radi- 
ation effects become more dominant; for strong enough 
magnetic field the particle loses most of its energy be- 
fore the next kicking. Then, the kick of the electric 
field "pump" energy into the particle, which again radi- 
ate most of this energy before the next kick. Obviously, 
if the kicking strength, K, is small, the particle will lose 
its energy in an exponential way, and will spiral in to 
zero velocity (in the velocity phase space). On the other 
hand, if K is strong enough, the energy gain due to the 
kicking is larger than energy loss due the radiation. In 
this case, the particle will increase its energy (not neces- 
sarily due to Arnol'd diffusion), and the velocity phase 
space will be limited. 

The interesting phenomena occur when the energy loss 
due to the radiation is approximately equal to the pump- 
ing energy due to the kicking of the electric field. In that 
case one can expect that : a) the velocity of the parti- 
cle will decrease to zero, b) the particle behavior (in the 
velocities phase space) is a chaotic behavior, c) the par- 
ticle will diffusion (while increasing its velocity) through 
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the web like structure, or d) the particle will increase 
its velocity in a stochastic way, not through a web like 
structure. 



B. The approximated solution for the strong 
magnetic field case 

For strong magnetic field, we assume T tq. In this 
case, it is more convenient to write Eq. ( p5[ ) as 

A) cosh(^- - Ci) = (1 + ^l-Pl) (1 + e) , (47) 



where 



Eq. 



£ = e T o 



can be written as, 



(48) 



a = Qt + Qtq In 



i + VT^ffi 



+ ftT ln(l + £ 2 ). (49) 



When T > To, ^ < 1 and in the first approximation, 
£ is negligible. Under those circumstances, Eqs. ( |26| ) 
become, 



1 + n 2 r 2 



n+1 



n+l 



(54) 



Note that since the velocities in Eq. ( p4| ) are limited 
by one, and since Oto <C 1 and To < 1 the difference 
x n+ \ — x n will be a very small number^. 

In order to observe some interesting behavior the ef- 
fective kicking should be strong enough. The increments 
in the phase are very small number since x n+ i — x n = 
0{tqVq) when B 3> 1. Thus, it is necessary to choose 
sufficiently large k in the kicking function, f{x n ) = 
K sin kx n (Eqs. ©-©), 



k = 



{—)■ 



In our calculations, we used 



TO WO 



(55) 



(56) 



Obviously, the behavior also depends on the kicking 
strength, K . For small kx the kicking approximated by 
f(x) = K sin kx w Kkx. Thus, in some cases, large K is 
equivalent to large k. 



dx _ 2 

~ i + vi - $ 



e T o x 



^] cos(fit + O ) 



dy = 2 
dt 1 + 



e T « x 



where 



— sm !!( 

#/ 



rir In 



cos(ftt + 0o ) 



1 + 



(50) 



(51) 



(52) 



In the present case (T ^> To) <fio <C 0, since VLT ~ 1 (see 
Eq. ©). 

Integration of Eq. (KG) yields, 



1 + ft 2 r 2 



m (i) + b 



'if' 1 " 



■xq. (53) 



VII. RESULTS - STRONG MAGNETIC FIELD 

In order to analyze the parameter space, one should 
first determine the most important parameters in the sys- 
tem. In our case the parameters are the kicking strength, 
K, and the magnetic field strength, B. The initial condi- 
tions, v x .o, Vyfi, xq, were chosen to be much smaller than 
the light velocity c, since in this regime the relativistic 
effects are negligible, and the radiation effects are iso- 
lated from the relativistic effects. However, in spite of 
those small initial conditions, one can easily determine 
what would happen if those initial conditions were in 
the order of the light velocity; for large enough kicking 
strength, K, the particle will accelerate to that regime. 

Generally, a parameter space mapping requires a cal- 
culation of Lyapunov exponents |l2| , |l3| . However, in our 
case, because of the decaying terms in Eqs. (|4l"|)-([f2"|), 
we are interesting in the decaying regions. For this, the 
standard deviation, 



Av — V< v 2 



> - < v 



(57) 



We now return to the initial charge e = — eo and thus 
chan ge f2 — + — Q and a — > —a. The discussions from Sec. 
IV C are valid also here; Eqs. (pl|)-p2]) remain the same 
and Eq. (43) should be replaced by 



'In order to get some impression about the size of the mag- 
netic field, consider the following values, to ~ 5 x 10 8 /B 2 , 
Qr w 9 x 10 15 /S, T/t ~ 7B x 1(T 16 . Since the term 
T/tq enters in exponential functions, T/tq > 3| is sufficient 
in our approximation. Thus the magnetic field should be 
B>5x 10 n T = 5 x 10 15 gauss. 



8 



is considered (after a sufficiently large number of iter- 
ations when the motion converges to its stable behav- 
ior). In the case of decaying velocities, v — > 0, and thus, 
Av — * 0. In the most cases of periodic motion v = const 
and thus, Av — 0. As will be shown in the following, 
there exists a bifurcation behavior in the system, and 
since the gap between bifurcation points decreases ex- 
ponentially, the period one region is much larger than 
the other periods (the exponential decrease start from 
the point between period one and period two, and do 
not include the period one itself). Moreover, since it is 
very difficult to balance the kicking strength, K, with 
the magnetic field strength, B, in such a way that the 
pumping energy due to the kick will be radiated before 
the next kick, just a small part of the parameter space 
will lead to a quasiperiodic or to a periodic motion. The 
dominant part of the parameter phase space should be 
either chaotic or stochastic, and thus should be indicated 
by a larger value of Av. A high value of Av simply means 
that the basin of attraction (if it exists) is larger. 

In Fig. ||a we present a density plot of the parameter 
phase space (the height is Av). The darker regions have 
higher Av. For small K the velocity decays to zero. For 
larger K the basin of attraction is larger, and one expect 
to find a more complex behavior. However, there exist 
stability lines which reflect that the particle either loses 
its energy or moves in a periodic motion. The stability 
lines seems to be continuous lines with well defined shape 
(the small discontinuities can be related to a high period 
level of the orbit; in this case Av but Av will be 
relatively small). Even for large K values, in which the 
particle is kicked very strongly, it can move in a periodic 
way. However, those lines accumulate just a minor part 
of the parameter phase space, and thus one can expect 
to find a chaotic behavior near the neighborhood of an 
arbitrary point in parameter phase space [ jl2| . Such a 
parameter phase space gives an indication of a robust 
chaotic system Jl| with three positive Lyapunov expo- 
nents!. Fi S- I 13 

shows a contour plot of Fig. km. 



§ According to the conjecture in Ref. | |12( |, if the number of 
positive Lyapunov exponents is larger than the number of the 
relevant parameters, the chaotic behavior is "robust"; for an 
arbitrary point in the parameter space one expects to find 
a chaotic behavior in the neighborhood of this point. When 
the number of positive Lyapunov exponents is less or equal 
to the number of relevant parameters, the chaotic behavior is 
less strong; in the neighborhood of an arbitrary point in the 
parameter space it is difficult to find a chaotic behavior. In 
our case, the map is a three dimensional map, and thus the 
fractal dimension must be less or equal to four (one should 
take into account the time, t; it did not enter explicitly in 
the map). In the more realistic case the map never covers the 
entire phase space and thus the fractal dimension will be less 
then four. Thus, the number of positive Lyapunov exponents 
should be three. 



* Itf* Fq. 5a 
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FIG. 5. The parameter space, B : K, mapping, (a) A 
density plot of the parameter space. The darker regions in- 
dicate high velocity standard deviation, (b) Contour plot of 
(a). 

Taking the a fixed magnetic field value, B, and chang- 
ing the kicking strength, K, will give more inset about 
the particle behavior. In Fig. ||a we have plotted a bi- 
furcation map using K as the relevant parameter. For 
each K value we have iterate the map for several thou- 
sand of iterations and plotted just the last few hundred 
of points (this have done to make sure that the particle 
was convergent to its basin of attraction). For small K 
values, the particle do not kicked strong enough, and it 
loses its energy exponentially to zero. For higher K val- 
ues (K = 1.35 x 10" 5 - 2.25 x 1(T 5 ) the particle has 
a constant velocity. Then, a series of bifurcation points 
appears, leading to the chaotic region (the enlargement 



9 



of the solid rectangle is shown in the inset, where the 
bifurcation points can be clearly seen). The chaotic re- 
gion which follows afterwards ends with other period one 
points, followed by another set of bifurcation points, and 
so on. The regular regions of Fig. ||a are parts of the 
stability line of Fig. ||. The maximum velocity that the 
particle reaches seems to grow linearly. 

Fig. 6a 
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Fig. 6b 
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FIG. 6. The bifurcation behavior of the radiative map (in 
the strong magnetic field limit), (a) The velocity v versus the 
kicking strength K. The inset shows an enlargement of the 
indicated rectangle, (b) The spacing between the consecutive 
bifurcation points, AK n , versus the bifurcation number, of 
two different cases (as mentioned in the text). The dashed 
line is the exponential fit of the two cases. 



Feigenbaum pfl . The bifurcation phenomenon is one of 
the signs of chaotic behavior. In Fig. ||b we present 
the exponential decaying behavior of the displacement 
between bifurcation points, AK ni in the inset of Fig. 
^a (case 1) as well as, for another bifurcation region 
(K = 8.25 x lCT 5 -8.45x 1CT 5 denoted as case 2). As can 
be seen, in the two cases the exponent is almost the same, 
B pa —1.545. The connection between the exponent B 
and universal constant S is given by, 



AK„ 



,1-545 



AK, 



n+l 



4.688 « S = 4.6692 (58) 



According to Feigenbaum O] the convergence to the uni- 
versal constant S should be when n — > oo (in AK n ). 
However, in the present case the convergence is rapid; 
the fitting of the exponential function is from the second 
bifurcation point. 

The evolution of the chaotic behavior as a function of 
the kicking strength, K, is shown in Fig. ^ (the initial 
conditions are V x $ = 0, V Vi q = 10~ 5 and the parameters 
values are B = 5 x 10 11 T, q = 4; each figure contains 
10000 points). In Fig. [7^, there are two separate regions 
which become closer for increasing K, and finally join to 
one region (Fig. |?]d). Moreover, the attractor occupies 
more area of the velocity space. The complexity of the 
attractor seems to be richer for larger K values. How- 
ever, the fractal dimension of the attractors of Fig. |t] do 
not increase monotonically (jl5| . The fractal dimension of 
Fig. [t] is summarized in Table |. We have calculated the 
fractal dimension, D2, using two different methods, the 
generalized information dimension (column 2 in Table |) 
and the generalized correlation dimension (column 3 in 
Table |) . The fact that the value of the fractal dimension 
is not an integer number is another indication for the 
chaotic behavior of the system. The higher the fractal 
dimension is, the greater is the complex behavior. Ad- 
ditionally, the non-integer fractal dimension indicates at 
least one positive Lyapunov exponent, (using the above 
arguments, there are three positive Lyapunov exponents) 
which is another indicator of chaotic behavior^ ■ 



*Tt follows from Kaplan- Yorke conjecture [|16| that if the 
fractal dimension is a non-integer number there exists, at 
least, one positive Lyapunov exponent. 



The series of bifurcation points is connected to the 
universal constant S ~ 4.6692, which was discovered by 
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Fig. 7 
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(d) 



-10 -5 5 10 -5 5 10 -5 5 

v,(x10" 6 ) 

FIG. 7. The chaotic attractors of increasing K values: 
A = 2.45 xl(T 5 , (b) A = 2.5 xl(T 5 , (c) K = 2.55 xlO -5 , 
if = 2.6x 1(T 5 , (e) K = 2.65 x HT 5 , (f) A = 2.7 x lCT 5 , (g) 
A' = 2.75 x 1(T 5 , (h) K = 2.8 x HT 5 , and (i) if = 2.85 x 1(T 5 . 
For larger A values the attractor occupied larger area of the 
velocity space. 

One of the signs of chaotic behavior is the folding- 
stretching phenomenon, as well as the self-similarity of 
the map 0. Fig. | shows the inner view of Fig. |b. 
Each figure is an enlargement of the previous indicated 
rectangle, and contains 10000 points. The self-similar 
structure is clearly seen from Fig. ||c and ||d. The loca- 
tion of the "attractor lines" seems to form a Cantor set 
which is usually observed in chaotic maps. The folding- 
stretching phenomenon which is connected to the self- 
similar structure can be seen, for example, in Fig. 
where the solid rectangle shows sharp turns of the maps. 
The Cantor set nature of the attractor lines is demon- 
strated in Fig. ^ where we choose a horizontal cut of 
Fig. ||a. There are more dense regions and empty re- 
gions, as expected. The inset of the figure shows that 
this behavior continues at smaller scales. 



Fig. # 


Info. D. 


Corr. D. 


a 


1.20 


1.13 


b 


1.32 


1.27 


c 


1.41 


1.39 


d 


1.43 


1.43 


e 


1.39 


1.41 


f 


1.36 


1.36 


g 


1.38 


1.32 


h 


1.47 


1.49 


i 


1.52 


1.53 



TABLE I. The fractal dimension, D 2 , of Fig. §. The first 
column is the figure number, the second column is the in- 
formation dimension, and the third column is the correlation 
dimension. 





-3.4 -3.2 -3 -2.8 -2.6 -2.4 
v x (x10" 6 ) 





1 











-2.88 -2.86 -2.860 -2.858 -2.856 -2.854 

v x (x1Q~") v x (x1Cf e ) 
FIG. 8. The self-similarity and the folding-stretching phe- 
nomenon (starting from Fig. ^b). Each figure is an enlarge- 
ment of the indicated rectangle of the previous figure. 



Fig. 9 
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FIG. 9. The Cantor set like behavior of the chaotic map. 
The inset shows the self-similarity of the attractor, keeping 
the Cantor set behavior of the lines. 

For large kicking strength the particle increases (or de- 
creases) its velocity to a region which balances between 
the kicking and the radiation. In Fig. [l^ an example 
for the limited map is shown (the parameters values are: 
B = 5x 10 u , q = 4, K = 30, v x . = 0, and v yfi = 0.99). 
The clear boundaries are caused by the radiation which 
decreases the velocity significantly before the next kick. 
The kicking pushes the particle toward light velocity; the 
gap between light velocity and the attractor boundaries 
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is approximately the energy loss caused by the radiation. 
The self-similar behavior and the strength-folding phe- 
nomenon can be seen in the inset of Fig. [l0| 

Fig. 10 




-0.5 



FIG. 10. The chaotic behavior of the radiative map (large 
B) for large kicking strength. 



VIII. INTERMEDIATE MAGNETIC FIELD 

In the previous sections the analysis of the weak and 
the strong magnetic field was discussed. In this section 
we will discuss the intermediate magnetic field strength 
regime. Basically, a numerical integration of Eq. ( prjj ) 
should be performed. However, the approximation for 
the weak magnetic field (Eqs. (pd|)-(|43|) is valid up to 
B = 10 10 T and thus one can avoid the numerical 
integration. 

In Fig. [ll] some of the possibilities that can appear 
in the intermediate magnetic field regime are presented. 
For example, Fig. 



Hi shows a destroyed web like struc- 



ture, i.e., the particle started to diffuse through a web 
structure, but after a sufficiently large number of itera- 
tions it either decays toward zero velocity, or stops its 
acceleration and iterates in a quasiperiodic motion (sec 
the enlargement of the indicated rectangle in Fig. Oa, 
where the particle decays towards a periodic orbit). The 
particle starts near an unstable fixed point (Eq. (|l6|)), 
and since the magnetic field is relatively high its influ- 
ence on the particle is larger and it decays inside. For 
smaller magnetic field the particle diffuses through the 



^As it already pointed out 



5 x W 8 /B 2 , Qt 



9 x 



10 15 /B, T/r « IB x 1CT 16 . Thus, the requirement B < 
10 14 gauss = 10 10 T will satisfied the condition T/tq < 1. 



filaments of the web. In might be that, if the initial 
condition were chosen closer to the unstable fixed point, 
the particle would diffuse. A supported example for that 
claim can be seen in Fig. [il"| b where the particle starts 
from smaller velocity (half of the velocity of Fig. |TT|a) us- 
ing the same parameters values as in Fig. [lT| a (and thus 
the approximation of Eq. ( fl6| ) is more accurate) and it 
diffuses to the high velocity regime (approximately, to a 
ten times larger velocity). Thus, one can safely conclude 
that for intermediate magnetic field the influence of the 
unstable fixed point is somehow weaker, and it is more 
difficult to observe the web structure. Or, to state it in a 
different way, the width of the filament is narrower when 
the magnetic filed is stronger. 

Fig. 11a 
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Fig. 11b 
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FIG. 11. The intermediate magnetic field regime. The 
parameter values are B — 5 x 10 6 T, q = 4, K = 10~ 5 , and 
Vx.o = 0. Each figure consist 10 s points. The initial velocity, 
Vyfi, in (a) is 2 x 10 -5 while in (b) is 1 x 10~ 5 . 
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IX. SUMMARY 

In the present paper we have investigated the effect of 
radiation on the stochastic web. Under the restriction of 
small magnetic fields (B < O(IOOT)) an iterative map 
was constructed. Moreover, the effect of radiation on the 
stochastic web is very small because of the small magnetic 
field. Qualitatively, the non-radiative and the radiative 
cases have similar web structure. Despite of the naive 
expectation that the diffusion rate in the radiative case 
should be smaller than non-radiative case, one can find 
cases in which the opposite effect is observed. 

Although it seems the effect of radiation is not qual- 
itatively significant under laboratory conditions, it have 
a strong effect in the presence of a strong magnetic field. 
Such a magnetic field could occur near or in a neuron 
star (or other heavy stars) and it can cause a large radi- 
ation correction to the motion of the particle. An iter- 
ative map was constructed for that regime also. It was 
found that the web structure disappears and is replaced 
by rich chaotic behavior which is demonstrated by the 
non-integer fractal dimension, bifurcation, self-similarity, 
and folding and stretching. The chaotic attractor occu- 
pies a small part of the velocity space, and practically 
means that the particle can not accelerate to infinite en- 
ergy, as in the regular stochastic web. This chaotic be- 
havior seems to be a "robust" chaotic behavior (implying 
that one can find chaotic behavior using slightly different 
parameter values). 

For the intermediate magnetic field strength, the influ- 
ence of the unstable fixed point is smaller relative to a 
weaker magnetic field. 
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